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Abstract Agricultural intensification in South Asia has 
resulted in the expansion and intensification of surface 
irrigation over the twentieth century. The resulting changes 
to the surface energy balance could affect the temperature 
contrasts between the South Asian land surface and the 
equatorial Indian Ocean, potentially altering the South 
Asian Summer Monsoon (SASM) circulation. Prior studies 
have noted apparent declines in the monsoon intensity over 
the twentieth century and have focused on how altered 
surface energy balances impact the SASM rainfall distri- 
bution. Here, we use the coupled Goddard Institute for 
Space Studies ModelE-R general circulation model to 
investigate the impact of intensifying irrigation on the 
large-scale SASM circulation over the twentieth century, 
including how the effect of irrigation compares to the 
impact of increasing greenhouse gas (GHG) forcing. We 
force our simulations with time-varying, historical esti- 
mates of irrigation, both alone and with twentieth century 
GHGs and other forcings. In the irrigation only experiment, 
irrigation rates correlate strongly with lower and upper 
level temperature contrasts between the Indian sub-conti- 
nent and the Indian Ocean (Pearson’s r = —0.66 and 
r = —0.46, respectively), important quantities that control 
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the strength of the SASM circulation. When GHG forcing 
is included, these correlations strengthen: r = —0.72 and 
r = —0.47 for lower and upper level temperature contrasts, 
respectively. Under irrigated conditions, the mean SASM 
intensity in the model decreases only slightly and insig- 
nificantly. However, in the simulation with irrigation and 
GHG forcing, inter-annual variability of the SASM circu- 
lation decreases by ~40 %, consistent with trends in the 
reanalysis products. This suggests that the inclusion of 
irrigation may be necessary to accurately simulate the 
historical trends and variability of the SASM system over 
the last 50 years. These findings suggest that intensifying 
irrigation, in concert with increased GHG forcing, is 
capable of reducing the variability of the simulated SASM 
circulation and altering the regional moisture transport by 
limiting the surface warming and reducing land-sea tem- 
perature gradients. 

Keywords South Asian Summer Monsoon • Irrigation • 
Latent heating • Monsoon intensity 

1 Introduction: intensifying irrigation and the South 
Asian Summer Monsoon 

Beginning in the 1960s, South Asia was witness to the 
Green Revolution: a period of unprecedented enhancement 
in agricultural production (Krishna Kumar et al. 2004) that 
helped avert potential regional famine conditions (Hazel 
2009). Expansion and intensification of irrigation was a 
significant part of this effort: 40 % of India’s agricultural 
production comes from irrigated crops and 80-90 % of all 
water withdrawals in India are for irrigation purposes 
(Douglas et al. 2006; Mall et al. 2006a, b; Lai 2011). A 
significant amount is withdrawn from groundwater 
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resources and introduced into the land-atmosphere domain, 
depleting these sources in many parts of India by up to 
20 cm/year (Douglas et al. 2009). Adding large amounts of 
irrigation water can alter the surface energy balance 
between latent and sensible heating, which then impacts 
surface temperature, atmospheric water vapor, regional 
circulation patterns and rainfall (Douglas et al. 2006, 2009; 
Saeed et al. 2009; Niyogi et al. 2010). In fact, the resultant 
cooling from irrigation could mask surface warming that 
may otherwise be attributed to increasing greenhouse gas 
(GHG, see Table 1 for a list of all acronyms used herein 
and their expanded forms) emissions (Bonfils and Lobell 
2007; Kueppers et al. 2007; Lobell et al. 2008). 

This is of particular concern for the South Asian Sum- 
mer Monsoon (SASM), which is dependent on the differ- 
ential heating between the land surface and the sea to 
establish its circulation in the late spring/early summer 
(Wang 2006). Furthermore, the most heavily irrigated areas 
are located in northern India, a particularly dry region that 
is subject to large temperature increases prior to the 
Monsoon Season (Niyogi et al. 2010). The pre-monsoon 
heating of this region is critical to help establish the pres- 
sure gradients that drive the initial monsoon circulation 
(Webster and Yang 1992; Wang 2006 and references 
therein). Has the intensification of irrigation in this region 
altered the SASM circulation? 

Some studies do suggest that the SASM intensity has 
weakened over 1960-2001 interval, perhaps due to reduc- 
tion in the upper-level and lower-level land-sea tempera- 
ture contrasts (Sun et al. 2010). This time period 
experiences both increasing irrigation and GHG forcing. 
GHG forcing has also contributed to rising equatorial 
Indian Ocean sea surface temperatures (SSTs), which have 
warmed by 0.5°-l °C compared to the early twentieth 
century (Alory and Meyers 2009). It is critical to under- 
stand the impact of each respective forcing (i.e., irrigation 
and GHG levels) on the monsoonal circulation. 

As the southwesterly monsoon winds converge moisture 
onto the Indian sub-continent, Monsoon Season rainfall is 
often used to assess changes in the SASM. Empirical and 


Table 1 Acronyms used and their expanded forms 


Acronym 

Expanded form 

GHG 

Greenhouse gas 

SASM 

South Asian Summer Monsoon 

GCM 

General circulation model 

MI 

Monsoon intensity 

CDC 

NCEP climate diagnostics center derived analysis 

CV 

Coefficient of variation 

LHF 

Latent heat flux 

ABCs 

Atmospheric brown clouds 


regional-model studies of the Indian- subcontinent have shown 
that intensive irrigation is more impactful to monsoonal rain- 
fall totals and distribution than the change from natural veg- 
etation to agricultural landscapes (Douglas et al. 2006, 2009; 
Lee et al. 2008; Niyogi et al. 2010). Increased latent heating 
over heavily irrigated regions alters the convective available 
potential energy and surface temperature, which may affect 
the SASM land-sea temperature contrasts (Douglas et al. 
2006, 2009; Lee et al. 2008; Saeed et al. 2009). Changes in the 
surface energy balance were also a consistent finding among 
general circulation model (GCM) experiments (Puma and 
Cook 2010; Cook et al. 2011; Guimberteau et al. 2011). This 
had the effect of delaying the SASM onset and weakening the 
Monsoon Season rainfall, which is consistent with regional 
model findings (Guimberteau et al. 2011) and highlights the 
potential impact of irrigation intensity to initial monsoon 
development. GCM simulations showed that as GHG forcing 
increased, an irrigation-induced cooling over South Asia was 
maintained, as high irrigation rates are able to meet increased 
evaporative demand (Puma and Cook 2010; Cook et al. 201 1). 
This is consistent with the ‘masking’ of increased temperatures 
due to GHG forcing over South Asia, which could impact 
SASM development in the future (Sen Roy et al. 2007, 2010). 
This surface cooling is also comparable to the observed 
warming of equatorial Indian SSTs and these combined effects 
could weaken the SASM circulation and moisture flux. Fur- 
thermore, it has been suggested that irrigation effects must be 
included in order to reproduce components of the SASM and 
its response to a warmer climate (Saeed et al. 2009). 

These irrigation-induced changes in South Asian surface 
energy balance may impact the large-scale monsoonal circu- 
lation, by way of a weakened land-sea temperature gradient. 
Many of the previous studies focus on the impact of irrigation 
on surface temperature and monsoonal rainfall — important 
quantities for many a livelihood in South Asia. However, an 
analysis of how intensifying South Asian irrigation have 
impacted the evolution of the large-scale SASM circulation 
over the twentieth century has not yet been attempted. Here, 
we conduct new GCM experiments to investigate the impact 
of irrigation intensity on the historical evolution of the SASM 
circulation. We focus on two research questions: (1) How has 
increasing irrigation affected the strength and variability of the 
twentieth century SASM circulation? (2) How do the irriga- 
tion impacts compare to those of increased GHG forcing and 
warmer regional SSTs? 

2 Methodology 

2.1 Model description 

For our simulations, we use the Goddard Institute for Space 
Studies (GISS) ModelE-R, an atmospheric GCM coupled 
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to the Russell-Schmidt dynamic ocean model (hereafter 
referred to as ModelE). ModelE was run at 2° x 2.5° 
(latitude x longitude) resolution with 40 vertical layers in 
the atmosphere and 32 layers in the ocean. The version of 
ModelE utilized here was frozen for IPCC Fifth Assess- 
ment Report, and builds upon the version of ModelE dis- 
cussed in detail in Schmidt et al. (2006) and Hansen et al. 
(2007), related to climate simulations from 1880 to 2000. 
This version also exhibits improvements in cloud param- 
eterizations and in the model code for sea-ice, ocean 
dynamics and the land-surface (Schmidt et al., in prep). 
The land surface is divided into vegetated and bare soil 
columns, extending to a depth of 3.5 m. Vegetation and 
phenology is prescribed (Matthews 1983, 1984) for eight 
vegetation types, distinguishing for its photosynthesis and 
stomatal conductance calculations using functions detailed 
in Farquhar et al. (1980) and Ball et al. (1987). Hansen 
et al. (2007) demonstrated that ModelE is able to reproduce 
the twentieth century climate, including trends, low and 
high frequency variability, when forced with observed 
SSTs. 

Wisser et al. (2010) estimated irrigation rates for the 
twentieth century as part of a reconstruction of twentieth 
century global hydrography. These authors produced their 
estimates by first using the University of Frankfurt/Food 
and Agriculture Organization Global Map of Irrigated 
Areas (Siebert et al. 2005a, b) to identify the areas equip- 
ped for irrigation at a spatial resolution of 5 arc minutes for 
the turn of the twenty first century (around the year 2000). 
From this map, they created a time series of irrigated areas 
in each grid cell extending back to 1901 by rescaling the 
year 2000 values with time series of irrigated area at the 
country level (compiled by Freydank and Siebert 2008). 
Prior to 1950, country level data on areas equipped for 
irrigation were unavailable for many locations. In these 
cases, Freydank and Siebert (2008) linearly extrapolated 
back to 1900. Using these irrigated area data, estimates of 
monthly irrigation rates were obtained using a water bal- 
ance and transport model (WBMplus) (e.g., Federer et al. 
2003). See Wisser et al. (2010) for further details. 

Water for irrigation is added to the top of the soil col- 
umn in the vegetated fraction of the grid cell (beneath the 
vegetation canopy). Surface water reservoirs (rivers and 


lakes) supply the initial irrigation in the same grid cell. In 
the event that the irrigation demand is not satisfied by these 
sources, water is then added to the system under the 
assumption that it is taken from groundwater sources that 
are not represented in the model. The flux of water is kept 
constant over the course of the day for days with non-zero 
irrigation, and is applied at every time step. Applied irri- 
gation water can either infiltrate the soil column (where it 
can then leave as evapotranspiration) or leave the grid cell 
as runoff. Additional details on the implementation of the 
irrigation scheme in ModelE can be found in Puma and 
Cook (2010) and Cook et al. (2011). 

2.2 Description of simulations and analysis 

We present results from four experiments performed with 
ModelE coupled to the Russell-Schmidt dynamic ocean 
model. Results from the four experiments are an ensemble 
average of five independent ensemble members all using 
the same respective forcings but with different initial 
conditions. For each member simulation, the model is 
started from an equilibrium state at 1850, following a spin- 
up period of 1,000 years. The simulations are analyzed for 
the 1960-2001 interval. The first simulation is a “Control” 
simulation in which all forcings are set to 1850 levels 
throughout the simulation and time- varying irrigation is not 
included. The second simulation, “GHG_only”, is per- 
formed with transient twentieth Century GHG and other 
forcings, but without time- varying irrigation. Also included 
are time- varying land surface specifications, however, we 
are not specifically assessing the effects of land-use 
change. The third simulation, “IRR_only”, includes time- 
varying irrigation over the twentieth Century, but all other 
forcings, including GHGs, are set to 1850 levels. Finally, 
the fourth simulation “GHG + IRR” applies both time- 
varying irrigation and all other transient twentieth Century 
forcings, including GHG concentrations and the resulting 
increases in equatorial Indian Ocean SSTs, which are 
comparable to the observed Hadley Sea Ice and Sea Sur- 
face Temperature (Shukla 2011). Table 2 provides a brief 
description of each model simulation discussed here along 
with their respective forcings. Figure 1 shows the ModelE 
specified irrigation requirements (or total irrigation water 


Table 2 ModelE experiments and description of forcings 

Simulation abbrev. 

Irrigation 

Forcings 

Control 
GHG_only 
IRR_only 
GHG + IRR 

No 

No 

Yes, time varying 
Yes, time varying 

GHG, solar, volcanoes and aerosols set to 1850 
Twentieth century GHG, solar, volcanoes and aerosols 
GHG, solar, volcanoes and aerosols set to 1850 
Twentieth century GHG, solar, volcanoes and aerosols 


All simulations underwent a spin-period of 1,000 years and were analyzed for the 1961-2001 time interval 
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specified in land- surface gridboxes from 5° to 40°N, 65 °- 
95 °E) over the Indian sub-continent during June-Septem- 
ber (Monsoon Season), and March- April-May (Pre-Mon- 
soon Dry Season). Totals are higher in the Monsoon 
Season, due to the greater availability of water via mon- 
soonal rainfall. The total irrigation is less in the Pre- 
Monsoon Dry Season, as most irrigation water is sourced 
from groundwater extractions, but the trend is increasing 
overall. 


3 Results 

3.1 The SASM circulation: ModelE representation 
and comparison of reanalysis products 

An important metric used in assessing SASM circulation 
strength is the Monsoon Intensity (MI) index (Sun et al. 
2010; Webster and Yang 1992). The MI index is defined as 
the difference between the average zonal winds at 850mb 
and 200mb over a region including the northern Indian 
Ocean and Indian Sub-continent (0-20N, 60E— 110E), 
computed for June-September. The strength of this shear- 
ing index is related to the intensity of the latitudinal tem- 
perature (and pressure) gradient between the monsoon 
trough region and the equatorial Indian Ocean. Figure 2a 
shows the 1960-2001 MI index as calculated from the 
ERA-40, NCEP twentieth Century and NCEP Climate 
Diagnostics Center (CDC) Derived Reanalysis products 
(Kalnay et al. 1996). A significant negative trend (assessed 
through a Mann-Kendall test at the 0.05 level) is found 
only in the NCEP CDC Derived product (Fig. 2b), while 
significant trends are not found in the ERA-40 or twentieth 
Century products. The NCEP twentieth Century reanalysis 
product is not as highly correlated with the previous two 


Irrigation over Time 



Fig. 1 Irrigation water (mm/day) added over the Indian sub-conti- 
nent (5°-30°N, 60°-95°E) in ModelE for Monsoon Season ( red) and 
pre-monsoon dry season {blue) between 1960 and 2001 


products. This discrepancy may be due to the exclusion of 
upper-level atmospheric temperatures in assimilation of the 
twentieth Century product (Compo et al. 2011), which are a 
critical component to the MI index calculation. 

ModelE reproduces the SASM’s characteristic south- 
westerly flow from the equatorial Indian Ocean and Ara- 
bian Sea during the Monsoon Season (Fig. 3). However, 
the model underestimates the strength of these winds. 
When forced with observed SSTs, the ModelE MI index is 
correlated to reanalysis-calculated MI indices, particularly 
the NCEP CDC Derived (Pearson’s r = 0.71). It is known 
that ENSO phases are related to a considerable amount of 
variability in the SASM as well, and the CDC Derived MI 
index has a 0.41 correlation with the SOI from 1960 to 
2001. ModelE (with prescribed SSTs) also produces a 
similar correlation to the SOI (R = 0.44). 

However, the coupled ModelE does not reproduce the 
full MI as depicted by any of the reanalysis products, and 
thus, produces an overall weaker monsoonal circulation 
(Fig. 2). This may be related, in part, to a cold bias that 
ModelE displays at higher altitudes over the northern 
Indian region relative to various reanalysis products, which 
can weaken the latitudinal temperature contrast that ini- 
tially drives the SASM circulation. Yet, the model is able 
to simulate the gross features of the SASM, and, most 
notably, it is able to simulate the magnitude and the overall 
variability of the All-India Rainfall (Parthasarathy et al. 
1995), shown in Fig. 4a. The model also reproduces known 
areas of maximum rainfall during the monsoon period, over 
the northern Bay of Bengal and over the Kerala coast 
(Fig. 4b, c). However, the model overestimates the rainfall 
in these areas and underestimates the rainfall in the interior 
of India. While the modeled total All-India rainfall is 
similar to the observed, the geographic distribution of 
rainfall is not fully reproduced. 

3.2 Monsoon intensity 

The ModelE simulations do not reproduce the full magni- 
tude or variability of the MI index shown in the various 
reanalysis products (Fig. 2a). When a Mann-Kendall test 
of significance was performed, no declining trend in MI 
was found in any of the ModelE simulations (Fig. 2b). 
Interestingly, prescribed (atmosphere-only) SST ModelE 
simulations with twentieth century GHG forcings and 
irrigation display significant declining trends (not shown 
here). While the GHG_only and Control simulations show 
an MI index with a greater standard deviation than any of 
the reanalysis products (see Table 3 for standard devia- 
tions), the GHG_only and Control simulations’ standard 
deviations were not found to be significantly different from 
the reanalysis products according to an F - test that checks 
for differences in variability. However, the IRR_only and 


5) Springer 


The response of the South Asian Summer Monsoon circulation 


25 


Fig. 2 a Monsoon intensity 
(MI) index (m/s) from 1960 to 
2001 (domain bounded by 
0-20°N, 60°-110°E) shown for 
the reanalysis products (ERA- 
40 — light blue , NCEP CDC 
Derived — red and NCEP 
twentieth Century — purple ) and 
ModelE simulations 
(GHG + IRR simulation — dark 
blue , GHG_only — green , 
IRR_only — orange and 
Control — gray), b MI index 
trend for 1960-2001 shown for 
each reanalysis product and 
ModelE simulation, in their 
respective colors. The 
denotes significant trends, 
evaluated by a Mann-Kendall 
test at the 0.05 level. The 
standard deviation for each 
series is shown in Table 3 


Monsoon Intensity Index 
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Trends in the Monsoon Intensity Index 
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denotes significant trend 


the GHG + IRR simulations did display reduced MI var- 
iability from the Control and GHG_only runs. 

To assess how MI variability changes with time, 10-year 
moving coefficients of variation (CV) are computed for 
each simulation’s MI index (Fig. 5a). The non-irrigated 
simulations distinctly differ from the reanalysis products 
and irrigated simulations, both in their range of CV and the 
trend in variability. The reanalysis products and ModelE 
simulations showed the following significant 1960-2001 
trends in their CVs (via a Mann-Kendall test at the 0.05 
level): CDC Derived = —0.056; twentieth Cen- 

tury = 0.119; ERA40 = -0.024; Control = -0.259; 
GHG_only = 0.097; IRR_only = 0.064; GHG + IRR = 
—0.135 (Fig. 5b). The irrigated simulations produce CV 


ranges more consistent with the reanalysis products. 
However, the GHG + IRR simulation, which also included 
the GHG warming equatorial Indian Ocean SSTs, is most 
consistent with all the reanalysis products in that it pro- 
duces both a significant declining trend and a more limited 
range in CV. 

3.3 Lower-level and upper-level land-sea temperature 
contrast 

The Monsoon Season lower-level and upper-level tempera- 
ture contrasts between the Indian sub-continent and the 
equatorial Indian Ocean are critical to the initial development 
of the large-scale S ASM circulation and MI. The upper-level 
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Fig. 3 850 mb Monsoon 
Season zonal winds (m/s) 
averaged over the period 
1960-2001, for a NCEP CDC 
Derived reanalysis and b the 
coupled ocean-atmosphere 
ModelE GHG_only run. In both 
maps, red areas (such as over 
the northern Indian Ocean) 
indicates westerly flow, while 
blue areas indicate easterly flow 


850mb Monsoon Season (JJAS) Zonal Winds (NCEP CDC Derived) 



850mb Monsoon Season (JJAS) Zonal Winds (GISS ModelE) 



temperature contrast is calculated as the difference in the 
200-500 mb thickness temperature between the northern 
Indian continent and equatorial Indian Ocean. The lower-level 
temperature contrast is calculated as the difference in the near 
surface temperatures over the northern Indian continent and 
the 500-800 mb thickness temperatures over the equatorial 
Indian Ocean (Sun et al. 2010). The northern Indian continent 
region is bound by 20°-40°N, 60°-100°E, while the equatorial 
Indian Ocean is bound by 10°S-10°N, 60°-100°E. Sun et al. 
(2010) found that the MI index was strongly correlated with 
the upper-level temperature contrast (Pearson’s r = 0.84) and 
moderately correlated with the lower-level temperature dif- 
ferentials (Pearson’s r = 0.51). The simulated MI indices for 
the GHG + IRR and GHG_only runs correlated with the 
upper-level temperature contrast by 0.76 and 0.74, respec- 
tively. For the IRR_only and Control runs, the correlations 


with the upper-level temperature contrast were 0.78 and 0.73, 
respectively. Correlations between the lower-level tempera- 
ture contrast and MI index were r = 0.53 for the Control run; 
0.27 for the IRR_only run; 0.34 for the GHG_only run and 
0.49 for the GHG + IRR run. Irrigation can change the sur- 
face energy balance, and thus are strongly (negatively) cor- 
related to the low-level temperatures. Correlations between 
the IRR_only and GHG + IRR irrigation and lower- and 
upper-level atmospheric temperature contrasts are noted in 
Table 4. The strongest correlation between irrigation and the 
lower-level temperature contrast exists during the Monsoon 
Season, although the correlation is also relatively high during 
the Pre-Monsoon Dry Season. 

The introduction of irrigation cools the land surface 
temperature in both the IRR_only and GHG + IRR sim- 
ulations when compared to the Control and GHG_only 
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All-India Monsoon Season Rainfall 




0 ° 30 ° 60 ° 90 ° 120 ° 150 ° 



< 


Monsoon Season Precipitation (mm/day) 


> 


10 


15 


20 


Fig. 4 a The All-India Monsoon Rainfall index (from 306 Indian 
weather stations, averaged over June-July-August-September) 
shown from 1960 to 2001 for the Merged CPC analysis (blue) and 
the ModelE GHG_only simulation (red). Monsoonal Season (June- 


July-August-September) Rainfall (mm/day) total averaged between 
1960 and 2001 for b merged CPC analysis and c the coupled ocean- 
atmosphere ModelE GHG_only 


simulations, respectively (Fig. 6a, b). The cooling is 
prominent over the (western) Indo-Gangetic Basin, corre- 
sponding to the most intensively irrigated areas. The 


300mb layer also displays a cooling (Fig. 6c, d) relative to 
the non-irrigated simulations, despite lesser correlations 
between the upper-level temperature contrast and surface 
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Table 3 Standard deviations of monsoon intensity (MI) index 


Time series of monsoon intensity index (m/s) 

Standard deviation 

NCEP CDC derived 

1.54 

NCEP twentieth century 

1.09 

ECMWF ERA-40 

1.44 

GHG + IRR 

0.89 

GHG_only 

1.69 

IRR_only 

0.92 

Control 

1.76 


irrigation. The cooling at altitude is more pronounced in 
the GHG + IRR simulation, which also includes the GHG- 
induced warming of equatorial Indian SSTs that raise 
atmospheric temperatures above the sea surface. Figure 7a 
shows the lower-level temperature contrast over time for 
each of the ModelE simulations. The IRR_only and 
GHG + IRR simulations display significant declining 
trends (assessed by Mann-Kendall tests at 0.05 level) with 
time (Fig. 7b), while no significant trends are shown in the 
non-irrigated simulations. Figure 7c indicates that as irri- 
gation increases, the lower-level temperature contrast 


Moving 10-Year Coefficient of Variation of Monsoon Intensity 
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Fig. 5 a 10- Year moving coefficient of variation (CV) in the 
monsoon intensity index from 1960 to 2001 (domain bounded by 
0-20°N, 60°-110°E), shown for the reanalysis products (ERA-40 — 
light blue , NCEP CDC Derived — red and NCEP twentieth Century — 
purple) and ModelE simulations (GHG + IRR simulation — dark 


blue, GHG_only — green, IRR_only — orange and Control — gray). 
b 1960-2001 Trends in the 10-year moving CV shown for each 
reanalysis product and ModelE simulation, in their respective colors. 
The denotes significant trends, evaluated by a Mann-Kendall test 
at the 0.05 level 
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Table 4 Pearson’s r correlation coefficient between irrigation requirements and atmospheric temperature contrast 


Simulation 

Pre-monsoon dry season 


Monsoon season 



Upper-level temperature 
contrast 

Lower-level temperature 
contrast 

Upper-level temperature 
contrast 

Lower-level temperature 
contrast 

IRR_only irrigation 
requirements 

-0.34 

-0.57 

-0.46 

-0.66 

GHG + IRR irrigation 
requirements 

-0.42 

-0.64 

-0.46 

-0.72 


decreases in both irrigated simulations. The lower-level 
temperature contrast is generally weaker in the 
GHG + IRR simulation than in the IRR_only simulation. 

Between the last and first decades of the GHG + IRR 
simulation, variability decreases by 37 % in the lower-level 
temperature contrast and 48 % in the upper-level temper- 
ature contrast. In the IRR_only simulation, the lower-level 
temperature contrast does not display a significant decrease 
in its standard deviation between the first and last decades. 
A 16 % decrease is produced between the first and last 
decades of the upper-level temperature contrasts. Only the 
GHG + IRR simulation displays significant decreasing 
trends (according to a Mann-Kendall test) in the 10-year 
moving CVs of the upper-level and lower-level tempera- 
ture contrasts (Fig. 8a, b), suggesting decreases in vari- 
ability since 1961. The decline is especially notable in the 
GHG + IRR lower-level temperature contrast. The 
IRR_only simulation, however, does not show the same 
degree of variation or significance of trend in the CV for 
the upper or lower-level temperature contrasts. 

3.4 Surface energy balance and moisture flux 

Large decreases in the Bowen ratio, representing increased 
latent heat flux (LHF), are shown over the most heavily 
irrigated areas in the IRR_only and GHG + IRR simulations 
(relative to their respective non-irrigated runs), particularly 
in northern India/Pakistan (Fig. 9a, b). The increase in 
modeled surface evaporation is comparable to the amount of 
irrigation water added (about 2 mm/day in May before 
SASM onset, for example) with negligible runoff (not 
shown). Such partitioning is consistent with the results for 
water-limited regions produced from studies regionally 
focused on South Asia, such Douglas et al. (2009). 

Increased LHF does not, however, necessarily result in 
an increased Monsoon Season rainfall. There is an increase 
in the simulated atmospheric water vapor over the most 
heavily irrigated areas in both irrigated simulations. The 
GHG + IRR and IRR_only simulations both show 
decreases in the meridional water vapor transport 
(Fig. 10a, b) over the northern Indian Ocean and Arabian 
Sea, which largely contribute to the Monsoon Season 
moisture convergence over the Indian continent. These 


decreases result in part from weakened lower-level wind 
fields in both simulations. 


4 Discussion and conclusions 

Agricultural intensification and increasing surface irriga- 
tion can be considered a forcing on regional climates. This 
study assessed the specific contribution of intensifying 
irrigation over the 20th century on the SASM circulation 
through analysis of the MI index and the driving land-sea 
temperature contrasts. Coupled GCM simulations were 
performed using time-varying irrigation specifications 
alone and in combination with GHG forcing. For each 
simulation, the magnitude and variability of the MI index, 
lower and upper-level temperature contrasts were evalu- 
ated along with the changes in the surface energy balance 
and moisture transport. The simulations were also com- 
pared with three reanalysis products, which were them- 
selves evaluated for mutual consistency. 

4.1 Impacts on the South Asian Summer Monsoon 
intensity 

ModelE is consistent with the ERA-40 and NCEP twentieth 
Century reanalysis products in that it does not produce a 
decline in the MI over the 1960-2001 interval. ModelE 
does not produce the full magnitude of the MI indicated in 
reanalysis products. However, the irrigated simulations 
produce MI variability that is more consistent with the 
reanalysis products and significantly decreased relative to 
the Control and GHG_only simulation. 

The variance in the MI index decreases at higher irri- 
gation rates (found towards the end of the twentieth cen- 
tury, as irrigation intensifies with time). The irrigated 
simulations show 10-year moving CVs that are more 
consistent with the three reanalysis products (all showing 
significant declining trends in variability), which suggests 
that the inclusion of irrigation may be necessary to accu- 
rately simulate the historical trends and variability of the 
SASM system over the last 50 years — a finding consistent 
with the regional climate studies performed by Saeed et al. 
(2009). 
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Fig. 6 Pre-Monsoon Dry 
Season (MAM) surface air 
temperatures (°C) for the 
a IRR_only and b GHG + IRR 
simulations, differenced with 
the Control and GHG_only 
simulations, respectively. Pre- 
Monsoon Dry Season (MAM) 
300mb temperatures (°C) are 
shown for c IRR_only and 
d GHG + IRR simulations, 
differenced with the Control and 
GHG_only simulations, 
respectively 
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The introduction of irrigation surface water dramatically 
alters the Bowen ratio. This change impacts the tempera- 
tures at the surface and at altitude over the most heavily 
irrigated regions, such as the Indo-Gangetic Basin. The 
build-up of the temperature and pressure gradients between 
this region and the equatorial Indian Ocean is critical to 
initiate and characterize the SASM circulation, particularly 
in the Pre-Monson Dry season. The irrigation-induced 
decreases in sensible heating limit the surface and altitude 
temperatures that initiate the low-pressure region that 
constitute the monsoonal trough. Irrigation effectively 
attenuates the SASM variability by inhibiting the rapid 
warming of the northern Indian sub-continent land surface. 
Although increasing irrigation impacted the modeled MI 
index variability, there is not a significant relationship 
between the absolute value of MI and the irrigation 


amount. Longer time series of the relevant SASM variables 
may be needed to assess the significance of these findings 
in the context of SASM variability. However, the impli- 
cations for less variable monsoons (particularly if rainfall is 
weakened) are relevant and extend to studies of predict- 
ability — a topic of prime importance in SASM research. 

4.2 Impacts to the lower- and upper-level land-sea 
temperature contrasts 

Significant trends were shown in the temperature contrasts 
between the Indian continent and the equatorial Indian 
Ocean, particularly in the lower-levels of the atmosphere 
and in the GHG + IRR simulation. As moderate-to- strong 
correlations exist between the MI index and the upper- and 
lower-level temperature contrasts, the variability in these 
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a Monsoon Season Lower Level Temperature Contrast 



Trends in the Lower Level Temperature Contrast 



C Lower Level Temperature Contrast vs Irrigation 



Fig. 7 a The Monsoon Season lower-level (800-500 mb) tempera- 
ture contrast (°C) between the northern Indian continent (20°-40°N, 
60°-100°E) and the equatorial Indian Ocean (10°S-10°N, 60°-100°E) 
from 1960 to 2001. Shown for the ModelE Control simulation (grey), 
the GHG_only simulation (green), the IRR_only simulation (orange) 
and the GHG + IRR simulation (dark blue), b 1960-2001 Trends in 


the lower-level temperature contrast shown for each reanalysis 
product and ModelE simulation, in their respective colors. The 
denotes significant trends, evaluated by a Mann-Kendall test at the 
0.05 level, c The Monsoon Season lower-level temperature contrast 
(°C) versus irrigation (mm/day) for the ModelE IRR_only (orange) 
and the GHG + IRR simulation (dark blue) 


quantities may also be related. The MI index’s decline in 
variability is consistent with the modeled trends in the 
temperature contrasts, particularly when GHG forcing is 


included. The lower-level temperature contrast, which is 
more directly impacted by changes in surface energy bal- 
ance, shows decreases in both its magnitude and variability 
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Fig. 8 Monsoon Season 
10-year moving coefficients of 
variation in the temperature 
contrast between the northern 
Indian continent (20°-40°N, 
60°-100°E) and the equatorial 
Indian Ocean (10°S-10°N, 60°- 
100°E) from 1960 to 2001 for 
the a lower-level (800-500 mb) 
temperature contrast and 
b upper-level (500-200 mb) 
temperature contrast. In both 
(a) and (b), the IRR_only 
simulation is shown in red, and 
the GHG + IRR simulation is 
shown in blue 


Fig. 9 Pre-Monsoon Season 
(MAM) Bowen Ratio for the 
a IRR_only simulation 
differenced with the Control 
simulation and b the 
GHG + IRR simulation 
differenced with the GHG_only 
simulation 
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Monsoon Season Meridional Water Vapor Transport 
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Fig. 10 Monsoon Season (JJAS) meridional water vapor transport ((m/s) * kg/kg) at 875 mb for the a IRR_only simulation differenced with the 
Control simulation and b the GHG + IRR simulation differenced with the GHG_only simulation 


under irrigation. Its decline may significantly contribute to 
limiting MI index variability. 

The yearly Monsoon Season LHF is strongly correlated 
to the lower-level temperature contrast (Pearson’s 
r = 0.85), highlighting the impact that changes in surface 
energy balance have on the monsoon circulation. The 
moderately strong correlations between the Monsoon 
Season lower-level temperature contrasts and the Pre- 
Monsoon Dry Season irrigation suggest that the dry-season 
irrigation partly inhibits the build-up of the lower-level 
temperature gradient. This may be particularly acute over 
the dry, water-limited northern Indian sub-continent, as 
water added for irrigation in the pre-monsoon season can 
evaporate efficiently, cooling the surface. Furthermore, 
although irrigation alone does not produce a significant 
trend in the upper-level temperatures, there is still a decline 
in the upper-level 10-year moving CV (though the impact 
is weaker compared to the low-level temperature con- 
trasts). The GHG + IRR simulation produced greater sig- 
nificant declines in the moving CV than IRR_only 
simulation, which indicates that GHG forcing, particularly 
in its warming of Indian Ocean SSTs, can act in concert 
with irrigation to further attenuate monsoon variability. 
Less variability and weakening of the temperature con- 
trasts, both at the surface and at altitude, are consistent with 
the shown decreases MI variability. 

4.3 Impacts to the atmospheric moisture transport 
and availability 

Under irrigated conditions, the major monsoonal circulatory 
components are made less variable when compared to the 


non-irrigated simulations. Both the IRR_only and 
GHG + IRR irrigated simulations produce weaker low- 
level southwesterly winds. This indicates that even irriga- 
tion alone can weaken components of the monsoonal cir- 
culation. Weaker SASM winds partially result from the 
significant decline in the lower-level temperature contrast 
and the cooling at upper atmospheric levels over the mon- 
soon trough region when intensified irrigation is present. 

Weaker SASM low-level winds would also impact low- 
level moisture transport, which is evidenced here. Atmo- 
spheric water content is increased over the heavily-irrigated 
Indo-Gangetic Basin. However, our results indicate that the 
southwesterly water vapor transport from the northern 
Indian Ocean and Arabian Sea is weakened under irrigated 
conditions, particularly in the GHG + IRR simulation. 
Coupled with reductions in atmospheric water vapor over 
the Indian peninsula, this result indicates that both increased 
irrigation and GHG forcing impact decreases moisture 
transport to the Indian sub-continent. Although slight 
declines in rainfall occur over the Kerala coast and northern 
Bay of Bengal, ModelE does not produce declines in rain- 
fall totals or variability that are significantly different from 
the observed All-India Rainfall average. These results are 
consistent with the findings of Douglas et al. (2009), and 
with the changes in the SASM circulation, specifically in 
the low-level southwesterly winds. 

4.4 Additional contributors to the South Asian Summer 
Monsoon variability 

This study demonstrated that the intensifying irrigation 
applications on the Indian sub-continent have potential to 
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limit the variability of the MI, decrease the magnitude and 
variability of the low-level driving monsoonal temperature 
gradients, and alter atmospheric water vapor transport and 
moisture fluxes. We also highlighted differences between 
the ERA-40, NCEP CDC Derived and the NCEP twentieth 
Century reanalysis products that should be considered 
when evaluating metrics such as the MI index. The NCEP 
CDC Derived product displayed a significant declining MI 
index trend that was not present in the other two reanalysis 
products. Furthermore, the NCEP twentieth Century 
product (which does not assimilate upper-level tempera- 
tures) did not produce the year-to-year variability that the 
other two products displayed 

It is important to note that the SASM system displays 
sensitivity to several other regional forcings and feedbacks 
that could further impact its variability under future climate 
conditions. Although we did not isolate and separately 
evaluate each of these forcings, atmospheric pollutants and 
particulate matter, fluctuations in snow accumulation on 
the Himalaya, land-use and land-cover changes (beyond 
irrigation intensification) and changes to natural drivers of 
regional climate variability (such as the El Nino Southern 
Oscillation) also contribute to the evolution and strength of 
the yearly monsoon. 

The South Asian region has recently been experiencing 
the effects of Atmospheric Brown Clouds (ABCs, air pol- 
lution clouds consisting of black carbon, sulphates, nitrates, 
fly ash, etc.) (Ramanathan et al. 2008). ABCs can warm the 
upper atmosphere, through high-altitude absorption of solar 
radiation, while cooling the surface. This can lead to fur- 
ther reductions in surface evaporation, shifting the energy 
balance and reducing the SASM latitudinal temperature 
gradient. Due to the widespread nature of ABCs at their 
peak, this change in energy partitioning can also decrease 
the vertical transport of water vapor and its availability 
over the ocean surface, which can lead to decreases in 
monsoonal rainfall (Cramer 2006). In contrast, other 
studies have indicated that high- altitude heating, particu- 
larly near the Tibetan Plateau, can intensify the heat low 
and associated pressure gradients responsible for driving 
the monsoon. Intensified monsoonal circulation may also 
draw a moisture flux off the ocean surface, which could 
lead to increases in monsoonal rainfall and changes to 
monsoonal variability (Lau et al. 2006). Additionally, 
stronger monsoons have been related to lower winter snow 
cover and higher winter minimum temperatures (Yanai and 
Wu 2006). As winter snow cover becomes more variable in 
a warmer climate, it may further increase SASM variabil- 
ity. The interaction between irrigation, ABCs and Hima- 
layan snowfall effects warrants further study to evaluate 
which is more influential in the region. 

In South Asia’s most heavily irrigated areas, ground- 
water and other external irrigation sources may cease to be 


economically viable in the early-mid twentyfirst century 
(Jha 2001). Decreased amounts of surface irrigation water, 
particularly during the dry season could allow for stronger 
surface warming, again increasing SASM variability and 
limiting water availability for existing rain-fed crops (as 
those reservoirs recharged by rainfall are now also used by 
farmers previously employing groundwater for irrigation). 

4.5 Summary of findings and statement of future work 

Previous studies have suggested decreases to the MI index 
over the 20th century. Irrigation decreased the simulated 
variability of both the MI index and the lower-level tem- 
perature contrasts, and it decreased the atmospheric mois- 
ture transport. When included, the GHG forcing amplified 
these decreases and further induced a declining trend in the 
variability of the driving temperature contrasts over the 
twentieth century. There were strong decreases to the 
Bowen ratio, which altered the surface energy balance, and 
transport of atmospheric water vapor. By including irri- 
gation along with GHG forcing, simulated MI variability 
was made more consistent with the reanalysis products. 
These findings are consistent with previous studies that 
evaluate the regional impact of irrigation on the moisture 
transport and rainfall associated with the SASM. This study 
further highlighted some inconsistencies between the 
reanalysis products used, which should be considered when 
they are used for future regional assessments in South Asia. 

Future work may focus on the relationships between 
altered surface energy balance, weakened low-level winds 
and regional moisture availability and transport. These will 
need to include the interactions with additional regional 
forcings, including those of “atmospheric brown clouds”, 
rising SSTs, changes to winter snowfall and rising surface 
temperatures outside of heavily irrigated regions. 
Groundwater used for extensive irrigation is also not an 
infinite resource, and future study building upon this work 
should account for changes in (or reductions to) irrigation 
water availability over the course of the twentyfirst century. 
Such studies will be useful to helping assess predictability 
of regional and local surface energy balance. 
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